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Abstract 

This  report  presents  numerical  procedures  for  solving  the  equations  of  Magneto-Fluid-Dynamics 
(MFD).  These  equations  consist  of  the  Navier-Stokes  equations  coupled  to  the  full  set  of 
Maxwell’s  equations.  They  govern  flow  within  electromagnetic  fields.  The  primary  interest  to 
the  Air  Force  of  this  work  is  flow  control  about  aerodynamic  vehicles  and  thrust  enhancement 
within  scramjet  engines  via  the  interaction  of  external  electromagnetic  fields,  produced  on  board, 
with  the  high  speed  flow  about  or  through  the  vehicle.  They  have  been  solved  for  both  external 
flow  and  internal  flow.  Test  cases  include  the  Hartmann  flow  problem^  weakly  ionized  flow 
simulating  the  RAM-C  flight  experiment^,  highly  ionized  flow  past  a  hemi-sphere-cylindrical 
body  simulating  the  Zieman  experiment^  and  flow  through  the  accelerator  section  of  and  “energy 
by-pass”  scramjet  engine"^’^. 


I.  Background 

The  interaction  of  an  ionized  flow,  which  caused  by  the  flow  crossing  the  bow  shock  of  a  high 
speed  aerodynamic  vehicle,  with  electromagnetic  fields  produced  on  board  offers  a  potential 
breakthrough  in  both  hypersonic  vehicle  design  and  propulsion.  Reductions  in  heat  transfer  and 
flow  control  using  magnetic  fields  can  be  important  for  enabling  a  hypersonic  vehicle  to  pass 
more  efficiently  and  safely  through  the  atmosphere.  Magnetic  and  electric  fields  placed  within 
the  propulsion  system  may  enable  the  extraction  of  electrical  energy  from  the  ionized  flow 
entering  the  engine,  while  simultaneously  slowing  the  fluid,  without  losses  in  total  pressure 
caused  by  shockwaves,  and  enhancing  complete  fuel  combustion.  The  extracted  energy  can  be 
returned  back  into  the  flow  after  combustion  for  further  flow  acceleration  and  engine  thrust. 
Realistic  aerodynamic  simulations,  under  the  conditions  of  expected  low  electrical  conductivities 
and  reasonably  strong  magnetic  fields,  will  be  required  to  determine  if  these  potential  benefits 
may  or  may  not  be  realizable.  Solutions  of  the  complete  equations  governing  magneto-fluid 
dynamics,  including  magnetic  induction  and  diffusion  within  strong  magnetic  fields,  are  needed 
to  perform  the  required  flow  simulations.  The  goal  of  this  research  has  been  to  develop 
algorithms  for  the  simulation  of  weakly  ionized  aerodynamic  flows,  both  internal  and  external, 
within  strong  externally  applied  electromagnetic  fields 

The  algorithms  developed  within  this  research  contain  numerical  procedures  for  solving  the 
governing  MFD  equations,  applying  boundary  conditions,  evaluating  temperature  and  pressures 
for  an  open  ended  set  of  chemical  species  composing  the  flow,  and  determining  the  electrical 
conductivity  of  the  ionized  gas.  The  governing  equations  of  MFD  consist  of  the  Navier-Stokes 
equations.  Maxwell’s  equations  and  the  equations  describing  chemical  and  thermal  equilibrium 
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and  non-equilibrium  high  temperature  air.  Maxwell’s  equations  appear  in  two  forms  -  1)  the 
usual  form  found  in  the  literature  and  2)  the  scalar  and  vector  potential  form  found  in  text  books, 
but  not  attempted  for  solution.  The  present  effort  developed  solution  and  boundary  condition 
procedures  for  both  forms  and  have  been  presented  at  seven  AIAA  (American  Institute  of 
Aeronautics  and  Astronautics)  meetings  during  the  last  three  years  '  .  This  volume  of  research 
cotributions  is  too  large  to  replicate  herein.  Instead,  the  procedures  will  be  briefly  described  and 
illustrative  applications  will  be  shown  from  the  AIAA  papers.  We  first  present  the  MFD 
equations. 


The  Equations  of  Magneto-Fluid  Dynamics 

XT  •  o  1  •  dU  dF  dG  dH  ^  .  .t 

1)  The  Navier-Stokes  equations - h - 1 - 1 - =  S  ,  with  U  =  \p,pu,pv,pw,e\  , 

dt  dx  dy  dz 

density  p ,  velocities  u,  v  and  w ,  total  energy  per  unit  volume  e . 

DE  1  r  V  X 

2)  Maxwell’s  equations  -  The  Ampere-Maxwell  equation  —  =  — - J 

St  Me  J 


Faraday’s  equation  —  =  -V x  E 

dt 


with  constraints 


V-E  =  —p^  and  V-5  =  0, 


and  charge  density  p" . 


electric  field  E ,  magnetic  field  B ,  current  density  J  ,  inductive  capacity  , 

magnetic  permeability  p^=4n:x  10"^ - -  and  charge  density  p" . 

{coulomb) 

3)  Generalized  Ohm’s  law  J  =  <j^{E  +  uxB) ,  with  electrical  conductivity  cr^ . 

4)  MFD  assumption:  charge  neutral  plasma,  p"  =  0. 


The  source  term  S  in  the  Navier-Stokes  equations  consists  of  the  Lorentz  force,  Lf  =JxB, 
acting  on  the  momentum  of  the  flow  and  Joule  heating,  caused  by  the  flow  of  electric  current 
through  the  fluid,  plus  the  magnetic  force  work  term  acting  on  the  energy  of  the  flow. 

S=  0,(JxB)^,iJxB)^,iJxB)^,-^J-J  +  iJxB)-u 

y  ^  G 

e 

The  Magnetic  Reynolds  number  is  defined  by  ,  where  Mq  and  /q  are  reference  flow 

speed.  Most  MFD  simulations  use  the  “Low  Magnetic  Reynolds”  number  approximation,  which 
assumes  that  the  electromagnetic  field  within  the  fluid  changes  insignificantly  from  that  imposed 
externally.  The  approach  taken  here  is  allow  the  fields  to  change  via  induction  as  prescribed  by 
the  full  set  of  Maxwell’s  equations.  It  will  be  shown  that  small  changes  to  the  electromagnetic 
field  can  produce  significant  changes  to  the  flow.  The  set  of  Maxwell’s  equations  above  are  the 
“usual”  set  found  in  the  literature.  This  set  consists  of  six  field  equations,  plus  the  constraint  that 

the  magnetic  field  remain  divergence  free,  V  •  5  =  0 .  The  “scalar  and  vector  potential”  form  of 
Maxwell’s  equations,  consisting  of  just  four  field  equations,  and  automatically  satisfying  the 
divergence  free  constraint,  is  given  below. 
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The  Scalar  and  Vector  Potential  Form  of  Maxwell ’s  Equations 

The  excellent  AIAA  paper  by  Giordano  describes  the  governing  equations  of  fluid  dynamics 
within  electromagnetic  fields.  His  paper  presents  and  derives  the  scalar  and  vector  potential  form 
of  Maxwell’s  equations,  summarized  below  following  his  description.  Because  the  divergence  of 

a  curl  is  zero,  it  is  desirable  to  represent  B  as  the  curl  of  a  vector  A,  thus  automatically 
satisfying  the  constraint  V  •  5  =  0 .  Substituting  B  =  Vx  A  into  Faraday’s  equation  yields 


=  -VxE  or  Vx  £  +  —  =0 


aVxA 


The  term  within  the  parenthesis  is  “curl  free’’,  which  implies  that  a  potential  (p  exists  such  that 

^  dA  -*  dE  ^  ^  1 

E-\ - =  -V (p .  Substituting  this  relation  into  Ampere’s  equation,  —  =  xB - J ,  and 

dt  dt 

using  the  identity,  VxB  =  VxVxA  =  V(V-A)  —  V^A,  yields 

d{V<p+^}  _  , 

dt  ^  ^ 


The  term  in  the  “curly’’  bracket  is  related  to  the  “Lorentz  gauge’’  and  is  set  to  zero  to  remove  the 

T  d 

arbitrariness  in  the  choice  of  the  vector  potential  A.  Thus  V-A  +  ^— =  0  and  the  time 

c?  dt 


d  A  22-  1  - 

dependent  equation  for  A  becomes  — —  =  c^V  A  A - j  .  From  the  constraint  on  the  electric 

dt 

-  -  1  -  -  SA  dV  •  A  1 

field,  V  E  =  —  p" ,  and  with  the  substitution  E  =  -Vg) - from  above,  - =  —  p'' , 

£e  dt  dt 

q2  ^2 

which  via  the  “Lorentz  gauge’’  becomes  — ^  =  c^V^cp  +  —  p’' . 

dt 

In  summary,  the  scalar  and  vector  potential  form  of  Maxwell’s  equations  are 


dt  8^ 


^  ^  ^  2^2  A  -*-  • 

cy^A+-j 

dt  £„ 


with  the  magnetic  field  B  and  electric  field  E  determined  from  the  solutions  for  scalar  (p  and 
vector  A  by 


B  =  S/xA 


E  =  -V  cp - 

dt 
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Note  that  only  four  elements  need  to  be  solved  for  in  time,  one  scalar  and  the  three  components 
of  A  and  also  the  constraint  V  •  5  =  0  is  automatically  satisfied. 

Solution  Procedures  for  the  Equations  of  Magneto-Fluid  Dynamics 

Decomposing  the  Electromagnetic  Fields  into  their  Imposed  and  Induced  Components 
It  is  advantageous  to  decompose  the  electromagnetic  fields  into  their  imposed  and  induced  parts. 
First,  the  induced  fields  may  be  orders  of  magnitude  smaller  than  the  imposed  fields  and  their 
magnitudes  may  therefore  be  compromised  in  numerical  precision  if  they  remain  combined 
during  numerical  operations.  Second,  the  implementation  of  initial  and  boundary  conditions 
becomes  simpler,  because  the  imposed  conditions  are  fixed  by  conditions  external  to,  and 
therefore  independent  of,  the  calculated  induced  flow  field.  Therefore,  we  write  for  the  magnetic, 
electric,  scalar  and  electric  potentials 

=B  +  Bq  and  =E  +  Eq 
(Pt=(p  +  %  and  4=A  +  ^ 

where  the  subscript  t  indicates  the  total  field  and  the  subscript  0  indicates  the  imposed  field.  The 
un-subscripted  variables  are  the  induced  field  components.  Note  that  both  forms,  the  usual  and 
the  scalar  and  vector  potential  forms,  of  Maxwell’s  equations  are  linear.  Therefore,  the  imposed 

fields,  Bq,  Eq,  (p^  and  ,  and  the  induced  fields,  B ,  E ,  (p  and  A ,  each  independently  satisfy 
Maxwell’s  equations. 

Boundary  Conditions  for  the  Electromagnetic  Eields 

Usually,  textbooks  give  electromagnetic  boundary  conditions  on  what  variables  or  derivatives 
are  continuous  across  boundaries,  so  as  relate  these  fields  inside  the  flow  field  to  those  outside. 
However,  it  is  desired  here  to  calculate  the  electromagnetic  fields  only  within  the  flow  field.  We 
therefore  need  the  values  of  the  electromagnetic  field  variables  or  their  derivatives,  the  so  called 
Dirichlet  and  Neumann  type  boundary  conditions,  instead  of  their  continuity  properties  across 
boundaries.  Otherwise,  both  the  region  within  the  flow  field  and  outside  would  need  to  be 
calculated.  To  avoid  this  complication,  we  assume  that  there  is  sufficient  information  concerning 
the  imposed  fields  along  all  boundaries  of  the  flow  field  to  set  the  needed  boundary  conditions. 
The  following  assumption  is  made  along  solid  wall  boundaries  for  the  applications  to  follow. 

Assumption:  We  assume  that  the  boundaries  at  solid  walls  are  “perfectly  conducting’’.  This 
implies  that  gradients  in  the  electric  field  components  parallel  to  the  walls  vanish.  Otherwise, 
infinite  currents  would  arise.  We  also  assume  that  these  components  of  the  imposed  electric  field 
are  constant  in  time  and  seek  solutions  converging  to  steady  state  in  time. 

Numerical  Method  for  the  Navier-Stokes  Equations 

The  Navier-Stokes  equations  contain  both  the  inviscid  Euler  and  viscous  terms.  All  terms  and  the 
boundary  conditions  are  treated  implicitly.  The  Euler  terms  use  a  modified  Steger-Warming  flux 
splitting  procedure.  The  approach  used  herein,  believed  originally  suggested  by  C.  Lombard^"^  in 
the  1980s,  has  the  benefit  of  partitioning  the  conservative  flux  vectors  themselves,  which,  for 
example,  unlike  the  state  vectors,  are  continuous  across  stationary  discontinuities  in  the  flow.  But 
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more  importantly,  this  approach  can  be  used  to  partition  general  flux  vectors  that  are  not 
homogeneous  of  degree  1  with  respect  to  the  state  vector  U  and  does  not  need  an  approximate 
Riemann  solution.  For  example,  the  flux  at  a  surface  between  mesh  points  i,  j  and  i  + 1,  j  is 
given  by 


F"  —  a" 


1+1/2, i 


F"+A"  F" 


(+1/2,; 


1+1,; 


where  JF  are  matrices  and  F",  and  F",  .  are  fluxes  at  and  about  the  flux  surface  i  +  l/2,  j . 

-(■+1/2,;  '’■/  '+'’•/ 

High  order  accurate  approximations  for  the  fluxes  can  be  made  by  upwind  extrapolation  or 
interpolation  of  the  flow  variables  to  the  flux  surfaces.  Second  and  third  order  accurate 
approximations  are  used  herein  along  with  a  TVD  (Total  Variation  Diminishing)  procedure. 


Numerical  Method  for  the  Maxwell  Equations 

(1)  The  Usual  Form  of  Maxwell’s  Equations  -  The  initial  attempt  to  solve  the  full  set  of 
Maxwell’s  equations,  using  a  similar  approach  to  that  for  the  Navier-Stokes  equations, 
showed  too  much  numerical  dissipation.  The  source  of  the  dissipation  was  found  and  a 
novel  scaling  of  the  Maxwell  flux  vector  was  able  to  remove  this  difficulty.  AIAA  Paper 

2008- 4010  describes  the  solution  procedure  for  solving  the  usual  form  of  Maxwell’s 
equations. 

(2)  The  Scalar  and  Vector  Potential  Form  of  Maxwell’s  Equations  -  Solution  of  this  form  of 
the  equations  was  attempted  because  of  the  “four  component”  form  of  the  unknowns  and 
the  automatic  satisfying  of  the  constraint  that  the  magnetic  field  be  divergence  free.  The 
four  component  form  is  fundamental  to  physical  laws,  i.e..  Momentum  and  energy,  where 
mass  is  a  form  of  energy.  Four  component  forms  transform  naturally  in  space-time,  i.e., 
under  Lorentz  transformations.  This  form  showed  no  solution  difficulties,  however 
boundary  conditions  were  more  complicated  This  work  is  described  in  AIAA  Paper 

2009- 0455 


Numerical  Method  for  the  Equations  of  Equilibrium  and  Non-Equilibrium 

An  air  chemistry  model  containing  nine  species,  plus  two  species  for  cesium  and  its  ion,  was 
considered.  The  cesium  species  were  used  to  seed  the  flow  within  a  scramjet  engine  to  enhance 
ionization  and  hence  the  electrical  conductivity  of  the  gas.  The  species  considered  were 
^2,02,  NO,  N,  O,  e~,  Nl,  NO\  N\  0\  Cs  and  Cs^ 


An  1 1  equation  reaction  model  for  the  1 1  air  species,  plus  cesium  and  its  ion,  is  shown  below. 


(1) 

N, 

+  M 

2N  +  M 

(2) 

0, 

+  M 

20  +  M 

(3) 

NO 

+  M 

N  +  O  +  M 

(4) 

N, 

+  0 

NO  +  N 

(5) 

NO 

+  0 

N  +  0^ 

(6) 

N 

+  0 

NO^+e- 

(7)  Cs  +  e  Cs  +  e  e 

(8)  N  +  e~  <r^N''+e~+e~ 

(9)  O  +  e~  +  e~  +  e~ 

(10)  V  +  V  V;+A 

(11)  0  +  6><^6>2+A 
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A  three  temperature  model  for  translational , rotational  and  vibrational  temperatures  was  used  to 
determine  the  thermal  non-equilibrium  of  the  gas.  This  work  is  deseribed  in  AIAA  Paper  2010- 
0225. 


Applications 

All  results  shown  below  were  obtained  with  the  numerieal  proeedures  diseussed  above.  The 
procedures  and  the  boundary  conditions  were  fully  implicit  and  run  to  high  CFL  numbers  (10^) 
limited  only  by  numerical  precision  considerations.  Each  problem  was  run  to  steady  state  in  less 
than  256  time  steps  on  a  lap  top  computer. 


1 
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Ram-C  Flight  Experiment  Figure  1  RAM-C  flow  field  pressure  contours 

During  the  1960’s  a  series  of  flight  experiments  were  made  at  high  altitudes  during  which 
electron  number  densities  were  measured  using  microwave  reflectometers.  The  vehicle  was  a 
sphere-cone  body,  with  a  9°  cone  half  angle  and  a  length  of  1.295m.  The  simulated  flow  field 
about  the  body  nose  is  shown  in  Figure  1  above,  showing  pressure  contours  for  a  free  stream 
(shown  in  solid  blue)  of  Mach  25.9  at  an  altitude  of  71km. 
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Figure  2  Temperature  contours,  RAM-C 


Figure  3  T  along  stagnation  streamline 


The  present  case  was  chosen  as  a  test  case  to  check  the  results  using  non-equilibrium  chemistry 
model  1  using  7  species,  N2,02,  NO,  NO^ ,N,0  and  e~ .  Cesium  was  included  but  played  no 

part  because  the  initial  mass  concentrations  were  set  to  zero.  A  single  energy  equation  was 
solved  for  temperature.  Temperature  contours  and  the  temperature  along  the  stagnation 
streamline  are  shown  in  Figures  2  and  3.  Species  mass  fractions,  =  pj  p ,  are  shown  in  Figure 

4.  Peak  electron  densities,  within  the  shock  layer  around  the  body,  are  compared  in  Figure  5.  The 
comparison  between  experiment,  shown  by  the  symbols,  and  computation  is  fairly  good.  This 
simulation  was  made  twenty  years  earlier  by  Candler^^  with  excellent  agreement.  Electron 
contours  are  shown  in  Figure  6.  Mass  fractions  contours  for  atomic  nitrogen  and  oxygen  are 
shown  in  Figs.  7  and  8. 


Figure  5  Peak  electron  number  density 
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Figure  7  Contours  of  N  mass  fractions 


Results  for  the  Ziemer  Experiment 


Figure  8  Contours  of  O  mass  fractions 
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Figure  9  Dipole  magnetic  field  lines 

In  1959  R.W.  Ziemer^  reported  results  from  an  experimental  investigation  in  magneto¬ 
aerodynamics.  He  placed  a  hemi-spherically  nosed  cylinder,  of  diameter  0.02m  and  made  of 
Pyrex  glass,  within  an  electromagnetic  shock  tube.  A  blast  wave  moved  at  Mach  21.5  into 
stationary  air,  at  temperature  273°  K  and  pressure  9.33N/m^,  past  the  model,  producing  a 
hypervelocity  flow  of  ionized  air.  The  period  of  steady  flow  time  was  only  about  10  to  20 
microseconds.  He  placed  a  coaxial  pulsed  magnetic  copper  coil  within  the  nose  of  the  body, 
which  produced  a  dipole  magnetic  field  of  4  Tesla  at  the  stagnation  point  of  the  body.  He 
observed  that  with  the  magnetic  field  turned  on  the  shock  wave  standoff  distance  increased  by  a 
factor  of  7.5  for  a  magnetic  interaction  parameter  2  =  ^ Poo^oo  =69.  The  free  stream 

conditions  for  the  present  Ziemer  flow  simulation  are  given  in  the  table  below.  Although  an 
attempt  was  made  to  match  the  experimental  conditions  of  Ziemer,  no  precise  match  for  the 
magnetic  interaction  parameter  Q  could  be  made  because  of  the  large  variation  in  electrical 
conductivity  cr^  behind  the  shock.  There  was  no  clear  choice  of  which  value  to  use.  The  imposed 

magnetic  field  at  the  stagnation  point  was  either  5o=0  or  5  Tesla.  The  temperature  behind  the 
shock  wave  near  the  nose  was  about  20,000°K  and  pressures  were  as  high  as  9x10"^  N/m^.  At 
these  conditions  the  flow  was  almost  completely  dissociated.  The  non-equilibrium  chemistry 
model  included  both  the  seven  and  eleven  species  model  for  non-equilibrium  air.  A  single 
temperature  was  used  in  the  present  simulation,  again  using  Park’s"^  thermodynamic  relations.  An 
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isothermal  wall  boundary  condition  was  assumed  with  the  wall  temperature  at  273°K.  However, 
a  catalytic  boundary  condition^^  for  both  species  and  internal  energy  was  used  at  the  wall. 


Free  Stream  Conditions 


velocity 

6715  m/s 

pressure 

5717  N/m^ 

temperature 

6049  °K 

density 

2.011xl0'^kg/m^ 

Figure  10  Pressure  contours,  Bo=0, 1  species  air 


Figure  11(a)  Pressure  contours,  5o=5,  7 
species  air 


Figure  11(b)  Pressure  contours,  5o=5,  11  species  air 


Standoff  Distance 

The  pressure  contours  for  the  flow  about  the  model  are  shown  above  in  Figs.  10,  11(a)  and  (b), 
with  and  without  the  dipole  magnetic  field  and  using  the  7  and  11  species  model  for  non- 
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equilibrium  air.  When  the  magnetic  field  was  turned  on,  the  standoff  distance  of  the  bow  shock 
wave  from  the  model  nose  increased  from  0.0017m  to  0.0058m,  a  factor  3.4  times  further,  for  the 
7  species  air  model,  and  to  0.0231m  for  the  11  species  model,  a  factor  7.7  times  greater  than  that 
shown  in  Fig.  10.  The  11  species  model  of  air  included  ionization  of  N  and  O  as  well  as 

NO ,  which  increased  the  number  of  free  electrons  and  therefore  the  electrical  conductivity  of 
the  gas  and  the  electromagnetic  field  interaction  with  the  flow.  The  11  species  air  model  is  more 
realistic  for  simulating  Zieman’s  experiment  because  the  temperatures  behind  the  bow  shock 
wave  reach  temperatures  as  high  as  20,000^  K. 

Imposed  and  Induced  Masnetic  Fields 

The  imposed  and  induced  x  components  of  the  magnetic  fields  surrounding  the  sphere-cylinder 
body  are  shown  below.  It  is  observed  that  the  induced  field  strengths  are  small  compared  to  those 
imposed  and  therefore  the  Low  Magnetic  Reynolds  number  approximation  would  probably 
suffice.  Also,  it  is  interesting  to  note  that  in  Figures  12  and  13(a)  and  (b)  the  induced  field  tends 
to  reduce  the  total  field  strength. 


0.03  p 
0.025  - 
0.02  - 
0.015  - 
0.01  - 
0.005  - 
>-  0  - 
-0.005  - 
-0.01  - 
-0.015  - 
-0.02  - 
-0.025  - 


-0.03 


U- 


J 


J 


-0.03  -0.02  -0.01 


0.01 


0.02 


0.03 


Figure  13(a)  Induced  field,  7  species  air 


Figure  13(b)  Induced  field,  1 1  species  air 
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Mass  Fractions  and  Temperature  along  Stamation  Streamline 

Mass  fractions  along  the  stagnation  streamline  are  shown  in  Figure  14  below  for  the  case  with 
the  magnetic  field  turned  on  for  the  seven  species  air  model.  The  plot  at  the  left  is  linear  and  that 
on  the  right  is  a  log  plot.  At  the  high  free  stream  temperature  diatomic  oxygen  is  near  completely 
dissociated. 


Figure  14  Mass  Fractions  pj  p  along  stagnation  stream  line,  left  -  linear  and  right  -  log  plots 


The  temperature  along  the  stagnation  streamline  is  shown  below  in  Figure  15,  again  for  the  7 
species  model  for  air. 


Figure  15  Temperature  along  stagnation  streamline 


Computational  Results  for  the  Hartmann  Problem 


Figure  16  Hartmann  Flow 


11 


An  appropriate  first  test  case  for  the  MFD  equations  is  Hartmann  Flow\  Hartmann  flow  consists 
of  viscous  flow  between  two  infinite  parallel  plates  within  an  imposed  constant  magnetic  field, 
normal  to  the  plates  and  an  imposed  electric  field,  E^<Q,  normal  to  the  flow  and  the 


dp 


imposed  magnetic  field.  A  stream  wise  pressure  gradient,  — ,  is  also  imposed  upon  the  flow.  All 

dx 

other  derivatives  with  respect  to  v ,  as  well  as  all  those  for  z  ,  vanish.  Analytic  solutions  exist  for 
velocity  u  and  induced  magnetic  field  component  as  functions  of  y  .  The  Hartmann  number 

is  defined  by  7?^  =  ,  where  //q  is  the  fluid  viscosity  and  y^  is  half  the  distance  between 


'Ao 


the  plates.  The  electrical  load  factor  is  defined  by  L  =  ■ 


The  mass  flow  rate 


r  yh 

udy  =  2u^yAs  held  constant  during  the  calculation,  i.e.  the  average  velocity  across  the 
channel  is  Mq.  Initially  the  current  has  only  a  z  component,  =<7^{^E^+uBy^ .  The  Lorentz 
force  acting  on  the  fluid,  L^=JxB,  has  x  component  Lj,^=—J^B^=—a^B^{^E^+uB^^^. 


Hence,  the  average  force  across  the  channel  is  =  -u^a^B 


- ^  +  1 

J 


.  Thus,  load  factors 


greater  than  one  will  accelerate  the  flow  and  those  less  than  one  will  decelerate  the  flow.  This 
force,  as  well  as  the  viscous  friction  forces,  must  be  balanced  by  the  pressure  gradient. 


dp 
dx 

viscous  forces. 


—  =  -u^crBl 
dx  ° 


-  +  - 


R. 


,  where  the  Hartmann  number  dependence  is  caused  by 


Computed  results  for  the  full  set  of  MFD  equations  are  compared  with  the  analytical  solution  in 
Figure  17.  The  initial  flow  was  at  Mach  0.1,  at  atmospheric  pressure  and  300“  K  temperature.  A 
perfect  gas  equation  of  state  was  assumed.  The  plates  were  isothermal,  also  at  300“  K,  and  the 
Hartmann  number  was  10.  The  electrical  conductivity  cr^  was  determined  from  the  Hartmann 

number  relation  given  above.  The  distribution  of  B^  is  anti-symmetric  across  the  channel  and  the 

velocity  u  is  symmetric.  The  numerical  solution  of  the  usual  form  of  Maxwell’s  equations, 
(filled  circle  symbols)  and  those  of  the  potential  form,  (open  square  symbols)  are  both  in 
excellent  agreement  with  theory*  (line).  The  results  shown  locked  in  after  128  time  steps,  but  the 
calculation  was  continued  to  2048  time  steps  without  further  change. 
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Figure  17  Hartmann  Flow  Solutions  Compared  with  Theory,  symbols-numerical,  curve-theory 


Computational  Results  for  an  MFD  Accelerator 


0.5 

y 

0 

- 

MFD  accelerator  section 

J _ 

-0.5 

1  1 

. 
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1  2  X 

3 

Figure  18  Sketch  of  MFD  channel  accelerator 

Park,  Bogdanoff  and  Mehta"^  presented  a  1-D  analysis  of  the  performance  of  a  scramjet 
propulsion  system  incorporating  the  MFD  (magneto-fluid  dynamics)  energy  bypass  concept^. 
The  system  contained  sections  for  an  MFD  generator,  combustor  and  an  MFD  accelerator.  The 
Park,  Bogdanoff  and  Mehta  accelerator  section  was  a  square  converging  duct,  2.846m  long,  of 
height/width  0.933m  at  the  entrance  and  0.730m  at  the  exit,  which  is  sketched  in  Figure  18.  It  is 
located  just  down  stream  of  the  combustor  section.  The  imposed  magnetic  field  across  the 
channel  was  =11 .28T  and  the  transverse  voltage  gradient  varied  from  £„  =  -30,990  V/m  at 

6  2 

the  entrance  to  -31,470  V/m  at  the  exit.  At  the  entrance  the  pressure  was  1.251x10  N/m  , 
temperature  3583°  K,  and  the  Mach  number  equaled  1.15.  The  electrical  conductivity  used  by 
Park,  Bogdanoff  and  Mehta  was  a^=  35. ?>lmho/m.  The  interaction  parameter 

2  =  =20  and  the  magnetic  Reynolds  number  =0.17,  defined  by  =  uj^cr^ju^ , 

where  Uq  and  1^  are  reference  flow  speed  and  accelerator  channel  length,  and  ju^  is  the 
magnetic  permeability. 


In  the  present  study,  the  flow  is  treated  as  two  dimensional  and  cr^  is  calculated.  Laplace’s 
equation  for  the  imposed  electro-magnetic  fields  within  the  channel  was  solved  to  keep 
VxBg  =  0 .  The  non-equilibrium  chemistry  model  contained  9  species, 
N2,02,  NO,  NO^ ,N,0  ,Cs,Cs^and  e~ ,  including  cesium  and  its  ion  with  initial  mass  fractions 

of  4x10'"^  and  0 ,  respectively.  Maxwell’s  equations  were  solved^  in  their  vector  potential  form. 
The  electrical  conductivity  depended  upon  the  degree  of  ionization  present  and  was  calculated 
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17  18 

via  Park’s  transport  model  .  The  rate  constants  for  cesium  came  from  Park  citing  the 
experimental  data  of  Chen,  C.J.,  Wu,  J.M.,  Wu,  F.T.,  and  Shaw,  D.T^^. 


The  components  of  the  induced  and  imposed  magnetic  and  electric  fields,  some  showing  field 
lines,  are  shown  in  Figures  19  to  24.  The  induced  fields  are  very  small  compared  with  the 
imposed  fields,  particularly  the  electric  field.  A  supersonic  flow  in  a  converging  channel  should 
slow  down,  but  the  results  shown  in  Figures  25  show  considerable  acceleration  of  the  flow, 
caused  by  the  interaction  of  the  electromagnetic  fields  with  the  ionized  flow.  These  results  are 
similar  to  those  obtained  for  flows  simulated  with  equilibrium  gas  chemistry^^.  The  maximum 
temperature  within  the  flow  was  about  4,400°  K  and  the  flow  reached  chemical  equilibrium.  The 
solution  converged  in  about  200  time  steps. 


Figure  19  Imposed  field 


Figure  20  Imposed  field 


Figure  21  Induced  B^  field 


0  1  2 
Figure  22  Induced  B  field 
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Figure  23  Imposed  field 


Figure  24  Induced  E  field 
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Figure  25  Velocity  vectors  and  u  velocity  contours 
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Thrust,  Velocity  and  the  Unstart  Problem 

Figure  26  shows  accelerator  section  thrust  and  maximum  velocity  along  the  channel. 


Thrust = 


[pu^+p] 

)ds 

[pu^+p] 

\ds 

and 


“max  =  — max  {M(x,y,z)} 


The  velocity  increases  by  a  factor  of  more  than  two,  in  agreement  with  the  load  factor  setting, 
and  the  thrust  by  more  than  50%.  However,  a  shock  wave  is  seen  moving  back  up  the  channel, 
located  at  about  x  =  -0.5  in  the  figure,  which  will  cause  the  engine  to  unstart.  The  design  Mach 
number  at  the  entrance  of  1.15  cut  it  too  close  and  a  larger  Mach  number  is  needed  as  a  safety 
factor  to  prevent  an  engine  unstart.  The  accelerator  section  of  the  channel  is  2.846m  long, 
starting  at  x=0.  Note  the  peak  in  velocity  at  the  end  of  the  section,  where  the  imposed  magnetic 
field  decreases,  as  shown  in  Figure  20.  Apparently,  there  is  also  a  defect  in  velocity  at  the  start  of 
the  accelerator  section,  where  the  imposed  magnetic  field  increases,  which  causes  the  flow  to  be 
subsonic,  starting  a  shock  wave  to  form,  which  moves  forward  as  observed  in  the  figure  above. 
This  phenomenon  was  also  observed  in  equilibrium  flow  simulations^.  The  explanation  for  this 
follows. 


Figure  26  Accelerator  Section  Thrust  and  Velocity 

The  Lorentz  force  acting  on  the  fluid  is  given  by  F^=JxB,  where  J  =  (T^(E  +  u  x B) .  If  E  =  0 

the  Lorentz  force  always  acts  to  decelerate  the  flow.  The  imposed  electric  field  is  chosen  to 
oppose  this  deceleration  and,  if  sufficiently  large,  accelerate  the  flow.  Examination  of  Figures  18 
and  21  shows  that  as  the  flow  approaches  the  accelerator  exit  at  x=2.846m  the  imposed  magnetic 
field  decreases  and  the  electric  field  remains  constant,  thus  favoring  acceleration.  After  passing 
through  the  exit,  the  electric  field  suddenly  vanishes,  thus  favoring  deceleration  and  explaining 
the  velocity  peak  at  the  accelerator  exit.  On  the  other  hand,  as  the  flow  approaches  the 
accelerator  entrance  at  x=0,  the  imposed  magnetic  field  increases  and  the  electric  field  is  zero, 
thus  favoring  deceleration.  After  passing  through  the  accelerator  entrance,  the  electric  field 
suddenly  increases,  thus  causing  acceleration  of  the  flow. 

Conclusion 

The  goal  of  the  present  research  is  to  develop  the  numerical  procedures  for  solving  the  equation 
of  magneto-fluid  dynamics.  This  work  is  important  to  the  Air  Force  to  determine  the  potential 
benefits  of  the  interaction  of  electromagnetic  fields  produced  on  board  an  aerodynamic  vehicle 
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with  the  high  speed  flow  about  the  vehicle  or  through  its  engines.  The  primary  application  of  this 
research  is  to  simulate  flows  necessary  to  design  hypersonic  vehicles.  The  equations  governing 
ionized  hypersonic  flow  include  the  Navier-Stokes  equations.  Maxwell’s  equations  and  the 
equations  of  chemical  and  thermal  non-equilibrium.  The  present  paper  summarizes  the  past  three 
year  effort  for  solving  these  equations.  The  details  of  this  research  are  contained  in  seven  papers 
presented  during  this  time  at  meetings  of  the  AIAA 

The  procedures  were  applied  to  the  RAM-C  flow  test  problem,  the  Ziemer  experiment,  the 
Hartmann  flow  problem  and  to  the  flow  within  an  MFD  channel  accelerator  containing  strong 
imposed  electric  and  magnetic  fields.  The  RAM-C  flow  case  demonstrated  the  accuracy  of  the 
non-equilibrium  flow  solver.  The  Ziemer  flow  simulation  agrees  with  the  experimental  finding 
that  the  standoff  distance  of  the  bow  shock  wave  can  increase  several  fold  because  of  the 
interaction  of  a  magnetic  field  with  the  ionized  flow.  The  Hartmann  flow  problem  verified  the 
accuracy  of  the  procedure  used  to  solve  Maxwell’s  equations. 

The  MFD  channel  accelerator  cases  demonstrated 

1)  the  capability  of  the  non-equilibrium  flow  procedure  to  include  cesium  as  a  seeding 
material  to  increase  the  ionization  of  the  air  gas  mixture  and  therefore  the  interaction  of  the 
flow  with  the  imposed  electromagnetic  field,  and 

2)  the  need  for  a  high  enough  supersonic  Mach  number  at  the  entrance  of  the  accelerator  to 
prevent  an  engine  unstart. 
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